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ABSTRACT 



Aims. Integral Field Spectroscopy (IFS) is a powerful approach for the study of nearby galaxies since it enables a detailed analysis of 
their resolved physical properties. Here we present the sample of nearby galaxies selected to exploit the two dimensional information 
provided by the IFS. 

' Methods. We observed a sample of 48 galaxies from the Local Universe with the PPAK Integral Field Spectroscopy unit (IFU), 

of the PMAS spectrograph, mounted at the 3.5 m telescope at Calar Alto Observatory (Almeria, Spain). Two different setups were 
used during these studies (low -V300- and medium -V600- resolution mode) covering a spectraf range of around 3700 - 7000 AA. 
(-H ' We developed a fufl automatic pipeline for the data reduction, that inciudes an analysis of the quality of the finai data products. We 

' applied a decoupling method to obtain the ionised gas and stellar content of these galaxies, and to derive the main physical properties 

I , of the gaiaxies. To asses the accuracy in the measurements of the different parameters, we performed a set of simulations to derive 

Q ■ the expected relative errors obtained with these data. In addition, we extracted two aperture, central and integrated spectra, from the 

^ ' datacubes. The main properties of the stellar populations and ionised gas of these galaxies and an estimate of their relative errors are 

[/3 I derived from those spectra, as well as from the whole datacubes. 

, Results. The comparison of the central spectrum extracted from the datacubes and the SDSS spectrum for those galaxies in common 

shows a good agreement between the derived values from both samples. We find differences in the properties of galaxies when 
comparing a central and an integrated spectra, showing the effects of the extracted aperture in the interpretation of the data. Finally, 
^ ^ we present two dimensional maps of some of the main properties derived with the decoupling procedure. 

^ ! 

CO ' Key words, surveys - techniques: spectroscopic - galaxies: abundances - stars; formation - galaxies: ISM - galaxies: stellar content 

oo ' 

1 . Introduction Spain). The wide field-of-view (FOV) of PPAK (74" x 65") al- 

lows us to sample galaxies in the Local Universe covering most 
O . Galaxies in the Local Universe are the final product of their t^eir optical extent. We have explored two different setups 
^ , cosmological evolution, and their actual properties are a direct (j^^ ^^^j medium spectral resolution) to finally obtain ~ 40000 
consequence of their assembling, star formation, enrichment and spectm, covering an area of ~ 50 arcmin^. To exploit at maxi- 
environmental effects (interactions/mergers). Since all these pro- ^um the capabilities of this insti-ument, we have observed galax- 
cesses are imprinted m the spectral features along their whole ex- j„ ^ jg^jghift range of 0.005 < z < 0.025, in a compro- 
tension, a two dimensional analysis can help us to understand the ^ise between the covered area and their sampling (spatial reso- 
H evolution that galaxies have undergone. The use of Integral Field i^tion) jhis redshift range corresponds to a luminosity distances 
Spectroscopy (IFS) enables us to build two dimensional maps of Dl = 19 - 106 Mpc, and one arcsec would correspond to a 
with a detailed characterisation of their spectroscopic proper- i^g^j g^^jg of ~ 92 - 527 pc, assuming a standard ACDM cos- 
ties, and therefore IFS is a powerful tool for the study of nearby oology ( Ho = 70.5, Q. = 0.27, A = 0.73. lHinshawitd]l2009h . 
galaxies. Thus, this study enables to distinguish structures in the galaxies 
In this paper, we present an exploratory study of 48 nearby (arms, bars, bulge/disk), although smaller interesting regions, as 
galaxies observed using the PMAS spectrograph (Postdam Multi individual H II regions, will not be resolved with these observa- 
Aperture Sp ectrograph, Roth et al. 2005) in the PPAK mode tions. 
jVerheiien et a l. 2004; Roth et al. 2005), mounted at the 3.5 m 

telescope of the Centro Astronomico Hispano Aleman (CAHA, Previous works have explored the use of different inte- 
gral field spectrophotometers for a detailed study of nearby 

* Based on observations collected at the Centro Astronomico galaxies. In particulai', the SAURON project (BaconetajJ 

Hispano Aleman (CAHA) at Calar Alto, operated jointly by the Max- 20011: Ide Zeeuw et al |2002|), and its extension Adas3D 



- '—I 

X 



Planck Institut fur Astronomic and the Institute de Astrofisica de (ICappellari et al.ll201 ll) are focused on the analysis of early-type 
Andalucfa (CSIC). galaxies and bulges of spirals at z < 0.01, to study their kinemat- 
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ics and stellar populations. Due to the distance of their objects 
(Dl < 42 Mpc) and the FOV covered by the SAURON insti-u- 
ment (33" x 41"), this study is mainly restricted to the central 
part of galaxies. The study of nearby spiral galaxies has been 
addressed by the VENGA survey (32 nearby spiral galaxies, 
Blanc et al.ll20Tol) using the VIRUS-P spectrograph jHill et al.1 
2008V the DiskMass Survey (Bershady et al. 2010a, 146 nearly 
face-o n galaxies) combining PPAK and the SparsePak spectro- 
graph ( lBershadvetal.ll2004i l2005h. and the PING S survey (17 
nearbv diskv galaxiesjRos ales-Ortega et al.ll201 0'). using PPAK, 
as in this work. In particular, the mosaicking designed for the 
PINGS survey allowed to map the H II regions along the whole 
extension their galaxies and to explore the two-dimensional 
metallicity structure of disks. Our sample comprises galaxies of 
different morphological types in the considered redshift range, 
with predominance of spiral galaxies. In addition, they were se- 
lected in size to fit in the FOV of the PPAK instrument, and there- 
fore, to map their physical properties in their whole extension in 
just one pointing. The results obtained from the PING survey 
can be compared with the work presented in this article to anal- 
yse the distance effects in determination and interpretation of 
their physical properties. 

The sample selection, observational strategy, data reduction 
and analysis of the accuracy of the final data set presented here 
formed part of the exploratory studies for the final selection for 
the optimal instrumental setup and sample of the CALIFA sur- 
vey ( Calar Alto Legacy Integra l Field spectroscopy Area sur- 
vexPl lSanchez et al.l201dl2011ah . CALIFA is an ongoing survey 
for observing ~ 600 nearby galaxies to explore the mechanisms 
who drives their evolution along the colour-magnitude diagram 
by studying their kinematics, ionised gas and stellar content. As 
part of this work, we have performed a wide set of simulations 
to check the quality of the data and the recovering of the phys- 
ical properties of the observed galaxies. Althoug this analysis 
has produced a bench mark for the observational strategy and 
data reduction scheme for future samples, e.g. in the CALIFA 
survey, the data described in this article comprise a well sepa- 
rated data set, which contains interesting physical information 
of the considered galaxies, and therefore we present them as a 
separated study. 

In this paper we focus on the sample (§|2]l, observations (§[3]), 
and data reduction (§|4]l. The first analysis and results from these 
data, including a detailed study of the accuracy of the derived 
physical parameters from the simulations, are presented in § |5] 
Further analysis for different aspects of galaxies, as the proper- 
ties of the ionised gas, and particular interesting objects will be 
presented in forthcoming papers. 



2. Sample 

Two different sample selections have been explored in this work: 
(1) galaxies extracted from the SDSS DR4 imaging sample 
brighter than r < 15.75 mag, at a redshift slice between 0.005 < 
z < 0.025 (selection in volume and limiting magnitude) and (2) 
face-on disk galaxies included in the DiskMass Survey (DMS, 
iBers hadv et al.l l2010al Verheijen, private communications) with 
appropriate sizes to fill the field-of-view (FOV) of the instru- 
ment (angular isophotal diameter selection). The first criterion 
is the standard method to select galaxies at a redshift range. In 
this case, the angular size of each galaxy is different and, there- 
fore the two dimensional information available for the sample is 
limited. On the other hand, the face-on galaxy sub-sample was 




Fig. 1. Subset of 9 of the observed galaxies, where the difference 
in projected sizes can be seen. All images are 1.5' x 1.5' (north 
is up, east to left). Images from SDSS survey. 



selected to solve this problem, but in this case the sample is not 
complete and representative of the Local Universe, and particu- 
lar care has to be taken when interpreting the results in terms of 
global evolution. 

In total, a subset of 48 galaxies were finally observed, ran- 
domly selected from those two samples depending on their vis- 
ibility at the considered nights. Four galaxies were partially ob- 
served and they are not included in this work. Table[T]lists the ob- 
served galaxies, including some basic information as the name, 
coordinates, a basic morphological classification (when avail- 
able), redshift, and V-band absolute magnitudes. Although there 
is a bias towards gas-rich late-type galaxies due to the selection 
of the sample, it is an heterogeneous sample that includes bulge- 
dominated and almost pure-disk galaxies, face-on and edge-on 
objects, and galaxies with and without bars. This is shown in 
Fig.[ri where some galaxies included in this study (images from 
the SDSS catalogue) are presented. 

Finally, we present in Figure |2] the colour-magnitude dia- 
gram of the observed objects, indicating those ones coming from 
each sub-sample. Regions coiTesponding to the red sequence, 
the g reen va lley and the_blue cloud are a lso indicated (e.g. , 
Stra teva et al.l 2001; B elletd.1 12004 l2006t iFaberetal.l 120071: 
Chen et al.ll20 10). Although the sample presented here is small, 
there is a significant number observed galaxies in both the red se- 
quence and the blue cloud, as well as the green valley. Therefore, 
this sample will allow a preliminary study of the global evolution 
of galaxies in the Local Universe. 



' |http : //www ■ caha . es/CALIFA/j 



3. Observations 

Observations were carried out during 15 nights in several 
observing runs. Information related to the weather condition 
during the observations are presented in Table |2l The full 
sample was observed at the 3.5 m telescope of the Calar 
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Galaxy 

(1) 


RA (J2000) 
(2) 


DEC (J2000) 
(3) 


Type 
(4) 


(5) 


R25 
(mag) 
(6) 


Mv 
(mag) 
(7) 


(B-V) 
(mag) 
(8) 


ObservingDate 
(yy-mm-dd) 
(9) 


Grid 
(10) 


Seeing 
(") 
(11) 


Airmass 
(mag) 
(12) 


KUG 1033+403 


10'' 36'" 12.4' 


+40°03' 12.6" 


S? 


0.0229 


0.56 


16.5 


0.87 


2009 


-03 


-23 


V300 




1.12 


2MASX J13062093+5318232 t 


13''06'"21.0' 


+53"I8'23.2" 


AGN 


0.0237 




15.9 


1.25 


2009 


-03 


-22 


V300 


2.0 


1.11 


2MASX J13193595+5330102 t 


13''19'"35.9' 


+53"30'09.8" 


IrS 


0.0248 




12.4 


0.72 


2009 


-03 


-24 


V300 


2.0 


1.13 


CGCG 045-001 


13''30'"43.5' 


+07"3l'48.8" 




0.0220 




15.4 


1.07 


2009 


-03 


-20 


V300 


1.3 


1.19 


CGCG 07 1-096 t 


13''00'"33.2' 


+ 10"07'47.8" 




0.0239 




15.1 


0.79 


2009 


-03 


-20 


V300 


1.1 


1.25 


CGCG 148-006 t 


07''44""57.4' 


+28°55'39.0" 




0.0234 




15.3 


0.75 


2009 


- 10 


-23 


V600 


1.1 


1.30 


CGCG 181-068 t 


09''37»M9.2' 


+33"49' 25.8" 




0.0227 




15.7 


1.17 


2009 


-03 


-20 


V300 


1.2 


1.06 


CGCG 213-041 t 


ll''21"M6.8' 


+40" 20' 43.4" 




0.0208 




15.5 


1.16 


2009 


-03 


-21 


V300 


1.7 


1.04 


CGCG 293-023 t 


02''30"'21.5' 


+56" 47' 29.5" 




0.0156 




15.6 


0.99 


2009 


-03 


-23 


V300 




1.18 


CGCG 428-059 t 


22" 14" 5 1.4' 


+ 13"42'54.9" 




0.0242 




15.7 


0.90 


2009 


-06 


-27 


V300 


1.0 


1.28 


CGCG 428-060 


22''14»'57.9' 


+ 13"45'42.9" 




0.0242 




13.7 


0.92 


2009 


-06 


-28 


V300 


1.2 


1.21 


CGCG 430-046 t 


23''00"M6.2' 


+ 13"37'07.9" 




0.0243 




15.3 


0.86 


2009 


- 10 


-23 


V600 


1.1 


1.10 


IC 2204 t 


07''41"18.F 


+34" 13' 55. 9" 


(R)SB(i)ab 


0.0155 


0.00 


14.6 


0.93 


2009 


- 10 


-22 


V600 


1.5 


1.12 


IC 2500 t 


09''42"23.4' 


+36"20'59.l" 




0.0221 




15.2 


1.17 


2009 


-03 


-21 


V300 


1.7 


1.14 


IC 2515 t 


09''54""39.4' 


+37"24'30.8" 


Sb 


0.0193 


0.72 


15.1 


1.11 


2009 


-03 


-22 


V300 


1.6 


1.05 


MRK 1477 t 


13" 16" 14.7' 


+41"29'40.l" 




0.0207 




15.4 


0.47 


2009 


-03 


-24 


V300 


1.8 


1.04 


NGC 99 t,* 


00"23"'59.5' 


+ 15"46'l3.7" 


Scd 


0.0177 


0.04 


14.0 


0.62 


2009 


- 10 


-22 


V600 


2.3 


1.11 


NGC 3820 t 


ll"42>»04.9' 


+ 10"23'03.3" 


S? 


0.0203 


0.17 


14.9 


0.95 


2009 


-03 


-22 


V300 


1.2 


1.40 


NGC 4109 t,* 


12"06"'5I.I' 


-i-42"59'44.3" 


Sa? 


0.0235 


0.05 


14.8 


0.98 


2009 


-03 


- 19 


V300 




1.05 


NGC 7570 t,* 


23" 16"' 44.7' 


+ 13"28'58.8" 


SBa 


0.0157 


0.23 


13.9 


0.81 


2009 


- 10 


-22 


V600 


2.0 


1.13 


UGC74 


00"08'°44.7' 


-i-04"36'45.l" 


SAB(is)c? 


0.0131 


0.05 


14.5 


0.85 


2009 


- 10 


- 17 


V600 


1.3 


1.23 


UGC 233 


00"24'»42.7' 


+ 14" 49' 28.8" 


SBbc D 


0.0176 


0.00 


14.6 


0.62 


2009 


- 10 


-23 


V600 


1.1 


1.20 


UGC 448 


00"42™22.0' 


+29"38'29.9" 


SABc 


0.0162 


0.04 


14.4 


1.16 


2009 


- 10 


- 19 


V600 


2.0 


1.08 


UGC 463 t 


00" 43" 32.4' 


+ 14"20'33.2" 


SAB(rs)c 


0.0148 


0.00 


13.5 


0.85 


2009 


- 10 


- 18 


V600 


1.8 


1.29 


UGC 1081 


0l"30"46.6' 


+21"26'25.5" 


SB(is)c 


0.0104 


0.03 


13.8 


0.83 


2009 


- 10 


-22 


V600 


2.4 


1.15 


UGC 1087 t 


0l"31'"26.6' 


+ 14"16'39.0" 


SA(rs)c 


0.0149 


0.03 


15.1 


0.69 


2009 


- 10 


- 17 


V600 


1.0 


1.13 


UGC 1529 


02"02'"31.0' 


+ ll"05'35.l" 


SA(rs)c 


0.0155 


0.07 


14.7 


1. 01 


2009 


- 10 


- 18 


V600 


1.5 


1.17 


UGC 1635 


02"08"27.7' 


+06"23'41.7" 


SAbc 


0.0115 


O.OI 


14.8 


0.96 


2009 


- 10 


-23 


V600 


1.2 


1.10 


UGC 1862 


02"24'"24.8' 


-02"09'44.5" 


SAB(rs)d pec? 


0.0046 


O.ll 


14.0 


0.83 


2009 


- 10 


-23 


V600 


1.3 


1.20 


UGC 3091 


04"33"56.l' 


+01"06'49.5" 


SAB(s)d? 


0.0184 




15.4 


0.88 


2009 


- 10 


- 17 


V600 


0.9 


1.40 


UGC 3140 


04"42'°54.9' 


+00"37'06.9" 


SA(i-s)c? 


0.0154 


0.06 


13.6 


0.97 


2009 


- 10 


- 19 


V600 


1.8 


1.36 


UGC 3701 


07"11'»42.7' 


+72" lO' 09.5" 


SA(rs)cd? 


0.0097 


0.00 


15.2 


0.87 


2009 


- 10 


- 17 


V600 


1.5 


1.30 


UGC 3997 ★ 


07"44'»38.7' 


+40"2l'58.5" 


Im? 


0.0197 


0.06 


16.0 


1.07 


2009 


- 10 


- 18 


V600 


1.5 


1.13 


UGC 4036 


07"51'»54.7' 


+73"00'56.5" 


SAB(r)b? 


0.0116 


0.06 


12.7 


0.98 


2009 


- 10 


- 19 


V600 


1.6 


1.38 


UGC 4107 


07"57"01.9' 


+49"34'02.5" 


SA(rs)c 


0.0117 


0.01 


13.9 


1.06 


2009 


- 10 


-22 


V600 


1.2 


1.13 


UGC 5100 t 


09"34"38.6' 


+05"50'29.9" 


SB(s)b 


0.0184 


0.22 


14.9 


1.02 


2009 


-03 


- 19 


V300 


1.1 


1.23 


UGC 6156 t 


ll"06'"31.3' 


+48"39'05.7" 


S? 


0.0245 


0.56 


15.1 


0.54 


2009 


-03 


-24 


V300 


1.5 


1.11 


UGC 6410 t 


11"24'"05.9' 


-i-45"48'39.9" 


SABc 


0.0187 


0.11 


14.8 


0.78 


2009 


-03 


-20 


V300 


1.5 


1.08 


UGC 7993 t 


12"50'"31.6' 


+52"07'22.5" 


Scd? 


0.0I6I 


0.90 


15.3 


0.79 


2009 


-03 


- 19 


V300 


1.3 


1.25 


UGC 9837 t 


15"23'"51.7' 


+58"03' 10.6" 


SAB(s)c 


0.0089 


0.02 


14.6 


0.57 


2007 


-06 


-23 


V300 




1.10 


UGC 9965 


15"40"'06.7' 


+20" 40' 50.2" 


SA(rs)c 


0.0151 


0.06 


14.4 


0.72 


2007 


-06 


-22 


V300 




1.39 


UGC 11318 


18"39'"12.2' 


+55"38'30.5" 


SB(rs)bc 


0.0196 


0.03 


14.1 


0.72 


2007 


-06 


-23 


V300 




1.10 


UGC 12250 t 


22"55'"35.9' 


+ 12"47'25.l" 


SBb 


0.0242 


0.22 


14.1 


1.04 


2009 


- 10 


- 18 


V600 


1.1 


1.33 


UGC 12391 


23"08"'57.2' 


+ 12"02'52.9" 


SAB(s)c 


0.0163 


0.04 


14.7 


0.77 


2009 


- 10 


- 17 


V600 


1.5 


1.20 



Table 1. General information about the observed galaxies. All the galaxies are included in 2MAS^ Galaxies marked with thave available 
SDSS spectroscopy. Galaxies marked with ★ have been observed with GALEjfl 

(1) Galaxy identification; Coordinates: (2) Right ascension and (3) Declination in J2000 Equinox from NEE|f|; (4) Morphological type from 
RC3 catalogu^ (5) Redshift from NED; (6) Rji- Radius where the surface brightness reaches yu^ = 25.0 mag/arcsec^, obtained from RC3 
catalogue; (7) My: Absolute B-band magnitude from NED; (8) (B-V); (9) Observing date; (10) PPAK mode utilised for each observation (see 
Table[3); (11) Average seeing value during the corresponding night observation. Data from the Calar Alto Observatory archive. (12) Airmass 
for the each observation. 



° Two Micron All Sky Survey (Skruts kie et al.ll2006h : http : //www . ipac . caltech . edu/2mass/ 1 
* The GALEX Ultraviolet Atlas of Nearby Galaxies fcil d e Paz et al.ii2007l) 
' NAS A/IPAC Extragalactic Database: |http:/ /nedwww . i pac ■ caltech . eduTI 
Third Reference Catalogue of Bright Galaxies: htt p: //vizier .u-strasbg ■ fr/viz-bin/VizieR?-source=VII/l 55 1 



Alto observatory with PMAS spe ctrograph (Postdam Multi 



Aperture Spectrograph. iRoth et a 



(IVerheiien et"ani2004t iRoth et al. 



1 l2005h in the PPAK mode 
2005V The PPAK fibre bun- 
dle consists of_382 fibres of 2.7" diameter each (see Fig. 5 
in 'Roth et al.l l2005l) . The science fibres (331 fibres) are con- 
centrated in a single hexagonal bundle, with a filling factor of 
~ 65 %. The sky background is sampled by 36 additional fibres, 
distributed in 6 bundles of 6 fibres each, distributed along a cir- 
cle ~ 90" from the centre of the instrument FOV. The sky-fibres 
are distributed among the science fibres within the pseudo-slit in 



order to have a good characterisation of the sky; the remaining 
15 fibres are used for calibration purposes. 

A dithering scheme with three pointings was adopted dur- 
ing the different observing runs, in order to cover the complete 
FOV of the bundle and to increase the spatial resolution of the 
data. This scheme has been already adopted in previous stud- 
ies using PPA K, with a considerable increase in the quality of 
the data (e.g., Sanchez et al."2007bVCastillo-Morales et al .l2010l: 
Perez-Gallego et al. 2010; Rosales-Ortega et al. 2010). The off- 
sets in RA-dec of the different pointings, with respect to the 
nominal coordinates of the targets, were (in arcsec): (1) (0, 0),(2) 
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-21 -20 -19 

SDSS r-band Abs. Mag. 

Fig. 2. u-r vs. M, colour-magnitude diagram of the observed 
sample. In blue circles we have indicated the DR4 selected sam- 
ple, and in green triangles those galaxies from the Disk Mass 
Survey (see text for more detail). The red sequence, green val- 
ley and blue cloud are noted with shaded colour areas. It is clear 
the observed sample still covers each region with a significant 
number of objects. 

Table 2. Atmospheric conditions during the observations. For 
each night, we present average values of the seeing, humid- 
ity, transparency of the atmosphere, and the extinction in the 
y-band. 



Observing 




Seeing Humidity 


Tranparency 


Extinction 


date 




(") (%) 


(%) in V— band (mag) 


2007-06-22* 


Clear 


56 


95 


0.22 


2007-06-23* 


Clear 


48 


95 


0.22 


2009-03-19 


Clear 


1.2 43 


95 


0.15 


2009-03-20 


Clear 


1.3 62 


85 


0.25 


2009-03-21 


Clouds 


1.7 71 


75 


0.28 


2009-03-22 


Clear 


2.0 70 


95 


0.15 


2009-03-23* 


Clouds 


47 


75 


0.25 


2009-03-24 


Clouds 


1.8 62 


65 


0.25 


2009-06-27 


Clear 


1.0 44 


95 


0.20 


2009-06-28 


Clear 


1.2 70 


90 


0.20 


2009-10-17 


Clouds 


1.2 79 


85 


0.25 


2009-10-18 


Clouds 


1.5 58 


70 


0.30 


2009-10-19 


Clouds 


1.8 75 


75 


0.20 


2009-10-22 


Clear 


1.6 86 


90 


0.20 


2009-10-23 


Clear 


1.2 50 


95 


0.15 


*Due to technical problems, there are not data from the seeing monitor for these nights. 


Table 3. PMAS/PPAK instrumental setups. 




Grid type 


V300 


V600 




Spectral coverage 


3620-7056 AA 


3845-7014 AA 




Dispersion 




3.2 A/px 


1.5 A/px 




FWHM" 




10.7 A 


5.4 A 




Exposure time * 


600 s X 3 


600 s X 3 





Notes. Spectral resolution (FWHM) for PPAK IFU with 150 fim 
fibers. **' Exposure time for each dithering pointing. 



(+1.56, +0.78) and (3) +(1.56,-0.78), which allows to cover 
the holes between fibers in the central bundle. The spatial re- 
composition of the three pointings is included in the standard 
data reduction scheme adopted for these data (described in §|4li. 

Two different setups were used during these observations, 
listed in Table [3] For the runs before October 2010, it was used 



the V300 grating, with a nominal spectral resolution of ~ 10.7 A 
(Full Width at Half Maximun -FWHM) and a wavelength range 
of AA3620 - 7056 AA. All the targets observed during this run 
coiTespond to the first sub-sample described in the previous sec- 
tion. After that date, the V600 grating was used, with a nomi- 
nal resolution of ~ 5.4 A (FWHM) and a wavelength range of 
/1/13845 - 7014 AA. Most of the targets observed during this run 
coiTespond to the DMS sub-sample. The main reason for that 
change was the upgrade of the PMAS CCD, from a 2k x 4k chip 
to a 4k X 4k chip, that nominally increases the wavelength range 
by a factor two at the same spectral resolution. This increase was 
hampered by a vignetting effect that affects to a ~ 30 % of the 
fibers/spectra, at the edges of the CCD (see Fig.|3]l. Fibers posi- 
tioning in the entrance slit are in such a way that the central and 
outer regions of the hexagon are located in the central regions of 
the slit, while the outer regions of the slit corresponds to an annu- 
lar ring. Due to this positioning of the fibers over the detector for 
this instrument, the vignetted regions correspond spatially to an 
annulus ring at about ~ 15" from the centre of the IFU. Despite 
this effect, all the fibers present unvignetting more than ~ 50 % 
of their spectral range. 

4. Data reduction 

The large amount of data generated in IFS observations made 
that the establishment of a precise and quick procedure for the 
data reduction was one of the key points to test during this 
study. The reduction was performed using a fully automatic 
pipeline, that operates without human intervention, producing 
both the scientific useful frames and a set of quality control mea- 
surements (on a web-based interface). The pipeline was based 
on the ro utines included in the R3D (Sanchez 2006) and E3D 
(ISanchez 2004) packages. The reduction consists of the standard 
steps for fibre-based integral-field spectroscopy. For the V300 
data, a master bias frame was created by averaging all the bias 
frames observed during the night and subtracted from the sci- 
ence frames. In the case of the V600, observed with the new 
PMAS CCD, which has basically no structure in the bias, a sin- 
gle value bias-level (derived on the basis of the over-scan val- 
ues), is subtracted to the frames. The different exposures taken 
at the same position on the sky were then combined, clipping 
the cosmic rays using IRAF0 routines. The next step was to de- 
termine the locations of the spectra on the CCD, using in this 
case a continuum illuminated exposure taken before the science 
exposures. Each spectrum was then extracted from the science 
frames. 

In order to reduce the effects of the cross-talk, we adopted 
a modified ver sion of the Gaussian-sup pression proposed by 
ISanchezI ('2006'), and fully described in ISanchez et all (1201 Ibh . 
This technique assumes a Gaussian profile for the projection of 
each fibre spectrum along the cross-dispersion axis, and it per- 
forms a Gaussian fitting to each of the fibres after subtracting 
the contribution of the adjacent fibres in an iterative process. 
The cross-talk is reduced to less than a 1 % when adopting this 
method. 

Each spectrum extracted along the disper sion direction i s 
stored in a row-stacked-spectrum file (RSS, ISanchezI |2004|) . 
Wavelength calibration was performed using HeHgCd lamp ex- 
posures obtained before and after each pointing, yielding an ac- 



- IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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4000 4500 5000 5500 6000 6500 7000 

Wavelength 

Fig. 3. Intensity map of the differential transmission fiber-to- 
fiber for the V600 grating. The drop of the transmission at the 
edges of the image is due to the vignetting. The ~ 15% of the 
CCD is affected this effect. 



curacy of ~ 10% of the nominal pixel scale (i.e., ~ 0.3 A for the 
V300 grating and ~ 0. 15 A for the V600 one). Differences in the 
relative fibre-to-fibre transmission throughput were corrected by 
comparing the wavelength-calibrated RSS science frames with 
the corresponding frames derived from sky exposures taken dur- 
ing the twilight. 

As mentioned before, PPAK is equipped with 36 fibres to 
sample the sky, distributed around the science fiber-bundle, in 
six small bundles of six fibers each one, at a distance of ~ 75 " 
from the centre of the FOV (seen Fig. 5 in lKelz et al.ll2006h . The 
objects selected to observe in these studies cover a substantial 
fraction of the FOV of the central PPAK bundle, but most of the 
sky fibers are free from emission by the targets. The procedure 
adopted to derive the night-sky spectrum was to combine the 
spectra corresponding to these fibers, performinga 2<j clipping 
rejection to remove any possible contaminatioifl Finally, the 
sky-spectrum was subtracted to all the spectra the correspond- 
ing frame. 

Spectrophotometric standard stars from lOkj(ll990h were ob- 
served during the nights to perform the flux calibration. The ex- 
tracted spectra of these stars were compared with the flux cal- 
ibrated spectra available in the webpage of the observatory^, 
obtaining the transformation function from observed counts to 
intensity. In this step, the routines in R3D package take into ac- 
count the airmass and extinction of the observation of both the 
spectrophotometric standard stars and the science targets. This 
procedure ensures a good relative flux calibration from the blue 
to the red part of the spectra, if the weather conditions along the 
night were stable. However, an absolute offset between the de- 
rived and real fluxes are expected due to the incomplete spatial 
coverage of the FOV in the observations of the calibration star 
despite the large fiber size of PPAK (2.7"), and second order ef- 
fects like small inaccuracies in the centering of the star and the 
effects of the variation of the seeing. 

After reducing the frames corresponding to each particular 
pointing, the science spectra corresponding to the three dithered 
exposures were combined in a single frame of 993 spectra. In 
order to take into account possible variations in the atmosphere 



^ We adopted this criterion since a 3cr clipping rejection does not 
remove the contamination by stars in the field, external regions of the 
galaxies, or a companion galaxy in IFU observations, due to the posi- 
tions of the fibers. 

* |http : //www ■ caha . es/pedraz/SSS/Oke/oke . html | 



transmission during the exposures, the procedure rescales the 
spectra to a common intensity by comparing the integrated spec- 
tra within an aperture of 20"/diameter, chosen in a compromise 
between the signal-to-noise and the depth of the images, and to 
avoid seeing effects. The final position tables for the individual 
fibers are generated taking into account their relative positions to 
the PPAK central bundle and the offsets provided for the dither 
scheme. The information included in the combined spectra and 
position tables are used to create the final datacube. 

Next step in the reduction process is to perform the dif- 
ferential atmospheric refraction (DAR) correction. To perform 
this coiTection, the pipeline creates an 2D image by co-adding 
the flux within 4500 - 5500AA (approximately the V-band), 
and looks for the peak emission within a region between 
[15:62,15:62] pixels. This centroid is used as initial reference, 
that will look for any shift in this centroid along the wavelength 
within a box of 6 x 6 pixels. Once the shift is determined, a 
polynomial order fitting of order 3 is applied to smooth the cor- 
rection, and the shift is applied wavelength to wavelength, re- 
centering the cube and correcting for the DAR. 

In order to get the best possible absolute flux calibration, 
we re-calibrated our data using SDSS photometry, whenever is 
available (37 targets, 33 of them with quality enough for the 
analysis). Of the five SDSS filters (ugriz), our spectrum cov- 
ers the passbands of two, namely the g (Agg - 4770 A) and r 
(/leff = 6231 A) filters. We measured the counts of each galaxy 
in these SDSS images inside a 30" diameter aperture. These 
counts were converted to flux following the counts-to-magnitude 
prescription in SDSS documentatiorQ We crosschecked this 
method measuring magnitudes of stars on each field and com- 
paring them with SDSS DR6 photometric catalogue, obtaining 
less than 0.05 magnitude dispersion for g < 17.5 mag. 

We then extracted the spectrophotometry from our reduced 
cubes, summing the flux of individual spectra inside a 30" di- 
ameter aperture around the peak intensity of the object, and con- 
volving this spectrum with the SDSS g and r filter passbands 
(Jester et al. 2005). Using these two data pairs, a scaling solu- 
tion of the form: 



FcM = Fo(A) X a 



(1) 



was derived for each band, where Fco„(A) and F(){A) are the re- 
calibrated and original fluxes, respectively, and a is the calibra- 
tion constant. The root-mean-square (rms) for these values is 
~ 4% in the g band and ~ 5%. Finally, the average of the two 
scaling constants derived for the g and r filters is adopted. In this 
way we were able to recalibrate the spectrum to match the SDSS 
photometry. 

Figure |4] shows two examples for the resulting re-calibrated 
spectra observed with the V300 and V600 grating. Together with 
the values for the g and r filters used in this re-calibration, we 
present the derived values for SDSS u and / bands. In addition to 
these measurements, the figure includes the fluxes corresponding 
to the Johnson B, V and Cousin Rc magnitudes determined on 
the basis of the SD SS photometry an d the transformation equa- 
tions presented by Jest er et aP (l2005h 0. There is a good agree- 
ment in all cases suggesting that this re-calibration method is 
reliable. 

Figure |5] (top panels) shows the g and r-band magnitudes 
derived from the datacubes and the SDSS images before (blue 



' http : //www ■ sdss . org/dr7/algorithms /f luxc al . html#counts2ma( 
' Summary of the transformations in 

|http : //www ■ sdss . org/ dr 7/ algori thms/sdssUBVRITransf orm . html | 
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4000 5000 6000 7000 

Wavelength [A] 




Wavelength [A] 

Fig. 4. Examples of the flux re-calibration based on SDSS pho- 
tometry for CGC071-96 (observed with the V300 grating, top 
panel) and CGCG430-046 (V600 grating, bottom panel). The 
values for the SDSS g and r filters used in the re-calibration are 
presented with blue solid circles, as well the derived values for 
SDSS u and i filters. Johnson B, V and Cousin Rq magnitudes 
are also shown (red solid squares) to test the flux calibration. 



solid circles) and after (red solid squares) the recalibration pro- 
cess. Both photometric data match for the 85% of the galaxies 
within a range of +0.15 mag before the recalibration (just 4 ob- 
jects present a higher dispersion below +0.25 mag). The possible 
differences in colours is explored in Fig. |5] bottom panel. Any 
possible error in the relative flux calibration from blue to red 
would produce a change in colour The dispersion (rms) around 
the one-to-one relation is ~ 0.1 mag, consistent with the ex- 
pected errors from propagating both the SDSS photometric er- 
rors (~ 4 - 5 %) and the absolute flux calibration accuracy of our 
data (~ 8 %). 

In the few cases (11 targets) when the SDSS photome- 
try is not available for a galaxy, we obtained broad-band im- 
ages using CAFOS at the 2.2 m telescope of the Calar Alto 
Observatory, Spain. The observations were caiTied out for 6 
nights from January 28th to March 30th, 2010, when images 
with the Johnson B, V, /, and Cousin Rq filters were obtained. 
The CCD camera was equipped with a SITe pixel chip. The pixel 
size was 24 yum, corresponding to 0.53". Sky flats were exposed 
in every dawn and dusk phase. Bias subtraction, flat-field divi- 
sion, and trimming were performed in the usual manner, using 
the IRAF package. 

We photo metrically calibrate d these im ages using GSC2.3 
dLaskeretal .12 008) and USNO-B (iMonet et al. 2003) catalogues 
as it follows. We obtained B and V photometry for our galax- 




0.4 0.6 0.8 1 

g-r SDSS 



Fig. 5. Top: Comparison of g (right panel) and r-band (left panel) 
magnitudes derived from the integrated spectrum within an aper- 
ture of 30"/diameter extracted from the V600 datacubes and the 
corresponding magnitudes derived from the SDSS images, be- 
fore (blue solid circles) and after (red solid squares) the flux cali- 
bration described in the text. Bottom: Comparison of g-r colour. 
Considering the eiTors in the absolute flux calibration (~ 8 %) 
and between g and r band (around 4-5 %), the expected eiTors 
in the colours are ~ 0.1 mag, consistent with the errors in SDSS 
photometric data. 



ies using the qphot routine of the IRAF and adjusted the zero- 
point (ZP) magnitude to the GSC2.3 and USNO-B photometry. 
For these ZP calibration we measured the magnitudes of sev- 
eral stars (7-8) in the same field of our galaxies and then we 
determined the difference between our magnitudes and the cata- 
logued ones. In general, the dispersion obtained for the ZP deter- 
mined using the USNO-B catalogue was nearly 0.4 mag on av- 
erage, so we decided to use the GSC2.3 ZP with a dispersion of 
0.2 mag on average. One ZP was calculated for each image with 
this method. Following the same procedure described before, we 
then convolved a similar aperture spectra extracted from the dat- 
acubes with the B and V filter responses, obtaining another pair 
of data sets for recalibrating these cubes. Considering that the 
stars from the GSC2.3 catalogue used in this work presented er- 
rors of 0.1 - 0.4 mag, we expect that our recalibration procedure 
ensures an absolute spectrophotometric calibration better than a 
~ 15 % in all the spectral range. 

Finally, the pipeline performs a set of quality control mea- 
surements, that are accessible to the user on a web-based inter- 
face, as indicated above. Among these measurements the most 
relevant are (i) cross-checks of the focus accuracy of the instru- 
ment, (ii) independent estimations of the accuracy of the wave- 
length calibration (by comparing the nominal wavelength of the 
most prominent emission lines included in the sampled wave- 
length range with the measured ones on the data, for each in- 
dividual spectra within the reduced RSS), (iii) estimations of 
accuracy of the sky-subtraction (from the residual flux of the 
strongest and more variable night-sky emission line in the con- 
sidered wavelength range -OI 5577 A), (iv) estimations of the 
night-sky brightness at the moment of the observation, (v) esti- 
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Table 4. Authomatic quality control measurements along the 
pipeline. The average 5cr limit V-band magnitude, the V-band 
magnitude of the sky and the accuracy of the wavelength caU- 
bration are presented for each night. 



Observing 


5tr limit 


Sky 


Wavelength 


date 


V-band 


V-band 


calibration 




(mag) 


magnitude 


(A) 


2007-06-22 


22.73 


20.50 


0.30 


2007-06-23 


22.50 


20.50 


0.30 


2009-03-19 


23.02 


20.25 


0.43 


2009-03-20 


23.02 


20.28 


0.34 


2009-03-21 


22.92 


20.10 


0.26 


2009-03-22 


22.84 


20.12 


0.39 


2009-03-23 


21.54 


19.07 


0.22 


2009-03-24 


23.19 


19.01 


0.32 


2009-06-27 


23.05 


20.13 


0.35 


2009-06-28 


23.75 


18.35 


0.30 


2009-10-17 


22.55 


20.05 


0.15 


2009-10-18 


21.58 


19.42 


0.20 


2009-10-19 


22.00 


19.83 


0.15 


2009-10-22 


22.42 


20.13 


0.15 


2009-10-23 


22.32 


19.73 


0.23 



mations of the depth of the final datacubes (by computing the 
signal-to-noise ratio in the spectral range 5670 - 5730 AA). For 
the observations presented in this paper, the average 5cr detec- 
tion limit is ~ 22.5 - 23.0 mag/arcsec^ (V-band) and the ac- 
curacy of the wavelength calibration is ~ 0.3 A. In Table |4] 
we present these values, together the V-band magnitude of the 
sky, for each observing night. On the basis of these quality con- 
trol procedures, it is possible to discard seven datacubes due to 
their poor quality. Most of these objects were observed under 
non photometric conditions, and the control has clearly identi- 
fied them to be discarded, as expected. 



5. Analysis and results 

In order to extract physical properties of the galaxy from the data 
set, it is necessary to identify the emission lines produced by 
ionised gas in the galaxies, and decouple this emission from the 
underlying stellar population. Particular care has to be taken in 
this decoupling technique, since some of the emission lines (e.g. 
Hj6) may be strongly affected by underlying absorption features. 
Several authors has developed different fitting codes to obtain 
information about the stellar content a nd star formation history 
of galaxies (e.g., Cid Fernandes et al. 2005; Ocvirk et al."2006'; 
ISarzi et all |2006; Koleva et al. 2009; MacArthur et al. 2009). 
In this paper we adop ted the decoupling method proposed by 
ISanchez et al.l (l201 Ibh . based in a full spectrum modelling. This 
method involves a linear combination of multiple stellar popu- 
lations and the non-linear effects of dust attenuation to obtain 
information about the ionised gas and the stellar content of the 
galaxies, and can be applied to each of the individual spectra of 
the datacubes to obtain the two dimensional distribution of the 
parameters derived with this procedure. 

In this section we present the method (§ 15. Il l and analy- 
sis of the accuracy of the physical properties derived with this 
method (§ 15.21 ). We present in § 15. 31 a comparison with the SDSS 
data, and we explore differences in the properties derived from 
two different aperture obtained from the data in § 15.41 Finally, 
a first analysis from the two dimensional spectroscopic data is 
presented (§ 15.51 ). 



5.1. Decoupling the stellar population and gas content 

The method for the decoupling of the stellar populations from 
the emission lines can be summarised as follows (for a more 
detailed description, see Sanchez et al. 2011b): (i) A set of de- 
tected emission lines was identified in the spectrum, (ii) The un- 
derlying stellar population was fitted by a linear combination 
of a sim ple grid of six s ingle- stellar populations (SSPs) tem- 
plates by VazdekisetalJ (|2010'), with three ages (0.089, 1.000 
and 17.7828 Gyrs) and two metalHcites (Z = 0.0004 and 0.03fl 
These templates were first coiTected for the appropriate systemic 
velocity and velocity dispersion (including the instrumental dis- 
persion), taking into account the dust attenuatior0. In addition, 
a spectral region of 30 A width around each detected emission 
line was masked prior to the li near fitting, including also the re- 
gions around bright sky-lines (ISanchez et al.ll2007ah . (iii) Once 
we derived a first approximation of the spectrum of the under- 
lying stellar population, this was subtracted from the original 
spectrum to obtain a pure emission-line spectrum, (iv) To de- 
rive the intensity of each detected emission line, this emission- 
lines spectrum was fitted to a single Gaussian function per emis- 
sion line plus a low order polynomial function for the local con- 
tinuum, (v) A pure gas-emission spectrum was created, based 
on the results of the last fitting procedure, using only the com- 
bination of Gaussian functions. The pure gas-emission model 
was then subtracted from the original spectrum to produce a 
spectrum of gas-free spectra, (vi) This spectrum is fitted again 
by a combination of SSPs, as described before (but without 
masking the spectral range around the emission lines, in this 
case), deriving the luminosity-weighted age, metallicity and dust 
content of the composite stellar population, (vii) Finally, indi- 
vidual emission-line fluxes are measured in each spectrum by 
considering spectral window regions around the most promi- 
nent emission lines (fitting of Gaussian funtions to [Oii]/l3727, 
U/3, [Om]/l/l4959,5007,HQ', [N ii]/l6548,83, [S ii]/l/l6717,6731). 
With this method, we obtain information about the ionised gas 
and the stellar content in the analysed spectrum. 

This decoupling method can be applied to any galaxy spec- 
trum, either integrated spectrum or individual spectra in the dat- 
acubes, to analyse the spatially resolved spectroscopic proper- 
ties. As an example, we show in Fig. |6] the best fitted gas- 
emission (red line) and stellar population (green line) spectra 
obtained when applying this decoupling method on the central 
spectrum of CGC07 1-096 galaxy (black line). The differences 
between the real spectrum and the final fitting are shown with a 
blue line. These residuals represent an eiTor of ~ 7% over the 
signal of the global spectrum, and it shows the good agreement 
of the global fitting. 

5.1.1. Star formation history 

To explore a more realistic interpretation of the computed stellar 
population, we decided to compare the galaxy data not only with 
SSPs, but also with spectra obtained from a given star formation 
history (S FH). With this aim, w e used the population synthesis 
models bv lVazdekis et alJ (l2010l) . These models provide the time 
(T) evolution of the synthesised spectra of galaxies characterised 
by an initial metallicity (Z), a star formation history (SFH) and 



' This grid of templates has been also adopted by IViironen et al.l 
( 201 1) in their analysis of UGC9837 galaxy from this data set. 

The dust attenuation is als o a free parameter i n this decoupling 
method. The extinction law of ICardelli et alj ( 119891) is adopted, with 
a ratio of total to selective attenuation Ry = 3.1 iJenldns|[l987il . 
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CGC071-096 




4000 4500 



5000 5500 
Wavelength (A) 



Fig. 6. Example of the decoupling of emission lines and under- 
lying stellar population for CGC07 1-096. We show the central 
spectrum in black line, the best-fit model for the emission lines 
in red, and the best fitted stellar population, modelled as a linear 
combination of different SSPs, plus the computed dust extintion, 
in green. Finally, the differences between the real spectrum and 
the final fitting is shown in blue. 



an initial mass function (IMF) {F ,i(T)\Z,S FHJMF 



sidered an pure exponential-declining SFH dSandagd 19861) : 



). We con- 



SFR (f) = exp(-(ro - 0/t) 



(2) 



where the steepness of the decay are regulated by a single pa- 
rameter T. To obtain the final spectra we considered a simple (not 
delayed) exponential, being Tq = 15.5 Gyr the age of the galaxy 
today, and a Salpeter IMF (a unimodal IMF with yU = 1.3). 
With these input data and using the tools provided in the MILES 
(Medium-resolution Isaac Newton Telescope library of empiri- 
cal spectra) websit^ we created a grid of spectra with solar and 
sub-solar metallicities (Z = 0.019, 0.0004) and t from 1 to 19 in 
steps of 2 Gyr, with a S/N ~ 50 per pixel in the continuum, and 
spectral resolutions corresponding to the V300 and V600 grism. 

In order to obtain the equivalence between the characteristic 
time scale parameter r and the stellar population parameters, i.e., 
luminosity-weighted age and metallicity, we applied the fitting 
procedure explained in the previous section, using as input the 
spectra generated from these SFHs. In Fig. |2]we represent the 
luminosity-weighted age and metallicity derived for the grid of 
spectra from the different SFH parametrised by r as a function 
of this parameter. As it can be seen, we obtain a clear dichotomy 
for the metallicity. This is expected since we are considering two 
different metallicities for the spectra created from a given SFH. 
In this sense, the metallicity is well recovered for all the values 
of T, although the metallicity is over-estimated in the case of the 
solar metallicity due mainly to the limited available values in the 
grid of models using in this experiment. 

In the case of the age, we obtain old stellar ages for low val- 
ues of T and younger ages until t > 10, w here the age becomes 
almost constant. This agrees with iGava zzi et al. (2002), who 
found that dwarf and normal elliptical galaxies, i.e. galaxies with 
older stellar populations, are well characterised by r < 3. They 
also found that galaxies with younger stellar populations present 
T > 4, as spiral galaxies (Sa-Sb) with values of 4 < r < 6, and 
blue compact dwarfs with t > 7. The equivalences obtained here 
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Fig. 7. Equivalence between the stellar population parameters, 
i.e the luminosity-weighted age (top panel) and the luminosity- 
weighted metallicity (bottom panel), and the characteristic time 
scale parameter t for the two SFH considered (red circle for solar 
metallicity, blue squares for sub-solar). We present the results of 
the fitting procedure applied to mock galaxies with the V600 
resolution, when the dust attenuation is fixed (filled symbols) 
and computed (open symbols). 



(Fig. |7]i will be used for a characterisation of the stellar popula- 
tions in §EM| 

5.2. Accuracy of the derived properties 

Follo wing the idea of recent works based on IFU observations, 
(e.g., iBershadv et ani2010bl: ISanchez et al.ll201 Ibh . we explore 
in this section the accuracy in the measurements of different pa- 
rameters of our interest (e.g., emission line fluxes, redshift, ve- 
locity dispersion, dust extinction, average stellar populations). 
With this aim, we created a set of mock spectra of galaxies with 
different kinematic galactic parameters (velocity dispersion, red- 
shift), underlying stellar populations (old or young SSP), ionised 
gas and signal-to-noise (S/N) in the continuum. These mock 
spectra were generated as follows. First, we consider an un- 
derlying stell ar population using the stellar population models 
provided by Vazdekis e t alJ (j2010|3 We consider four possibil- 
ities for the underlying stellar population: a single stellar pop- 
ulation (young -0.7 Gyr and solar metallicity (Z = 0.019) - 
or old -12.6 Gyr and metallicity of Z = 0.004), or two stellar 
populations mixed in different ratios (10% of young plus 90% 
of old stellar population -2SSP-10-, and 1% of young plus 99% 
of old stellar population -2SSP-01). Next, we add an emission- 
line spectrum to simulate the ionised gas in the mock galaxy. To 
obtain more realistic data, we measured the emission lines on 
several regions of NGC628 from the 2D spectroscopic data pre- 
sented bv iSanchez et al.l (1201 Ibl) . and we scaled these values to 
obtain similar values than the ones measured in a real galaxy of 



|http : //www ■ iac . es/proyecto/miles/pages/webtools/get- sp^t-pitittfQit'AflftBfflatplqsl/proyecto/miles/ 1 
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Fig. 8. Relative errors in the Hafli/3 ratio when the dust atten- 
uation is fitted (blue squares) or fixed (red circle) for different 
equivalent width of the H/? emission line (EWh/s), for a spec- 
trum of a mock galaxy with a pure old (left panel) or young 
underlying stellar population (as explained in the text). 



our sample (see e.g. IViironen et al.]|201 The nebular contin- 
uum is considered negligible. Then we apply different values of 
dust attenuation, same for the stellar continuum and the emission 
lines. Next, we consider the kinematics of the mock galaxy, i.e., 
we broad the spectrum to a given velocity dispersion cr and shift 
it to a given redshift z- Finally, we add noise to the spectrum to 
simulate a real galaxy observation. As result of this process, we 
obtain a spectrum of a mock galaxy with well-known parameters 
that we can change to check the recovering of the information for 
these studies. 

In next sections we present the different studies that we have 
carried out to analyse the accuracy in measurements of the dust 
attenuation (§ 15.2.11 ). kinematics (§ 15.2.21 ). ionised gas (§ 15.2.31 ). 
and underlying stellar population (§ 15.2.41 ). of the galaxies ob- 
served with the two setups (V300 and V600) used in this work. 

5.2.1. Dust attenuation 

The treatment of the dust attenuation may affect the re- 
sulting derived parameters of the stellar p opulation (i.e. , 
the luminosity-weighted age and metallicity, ICalzettil 120011: 
[Munoz-Mateos et al. 2007, 2009). As pa rt of the decouplin, 
method explained above, we adopted the iCardelli et al.l (Il98' 
law, which may be not the optimal sol ution to study t he dust at- 
tenuation in star-formin g galaxies (e.g.. lCaIzettil200lT) . Recently, 
iMacArthur et al.l 12009) considered the attenuation law based on 
the two component dust models by Chariot & Fall (2000), es- 
pecially developed for star-formation galaxies. Despite the con- 
ceptual differences between these attenuation laws, both of them 
present a very similar shape in the wavelength range of our data, 
and it is not wide enough to distinguish between these different 
laws. 

To explore the impact of fitting the dust attenuation on the 
derived parameters, we generated a set of mock spectra of galax- 
ies following the procedure explained above. Since values of 
equivalent widths of the H/3 emission line EWf^ > 3.5 A are 
suitable to perform an abundance analysis in galaxies despite the 



dependence on other physical properties of the ionising cluster 
(e.g., its age), we decided to explore a range of equivalent widths 
of EWh/s ~ 1.5 - 15 A for our mock galaxies. Then we applied 
the fitting technique (i) fixing the value for the dust attenuation, 
or (ii) computing Ay with the fitting program, with values of Ay 
varying between Ay = 0.1 and 1 .2, in steps of 0. 1 . 

We present the results of this experiment in Fig. [8] where 
we show the relative random eiTors obtained for the Ha/Hfi ratio 
(in absolute values) when the dust attenuation fixed (red circles) 
and fitted (blue squares), as a function of EWh/j. In all cases, 
the HafH/3 ratio is recovered with a relative error < 20 % for 
EWh/3 > 8 a for an old underlying stellar population (left panel), 
while the relative errors for a young stellar population is < 20 % 
for EWh/s > 3.5 A (right panel). Although there are not dif- 
ferences in the determination of HafH/3 for an old population 
depending on the method (fitting or not the dust attenuation), 
slightly higher eiTors are found when the dust attenuation is fitted 
in the case of a young underlying stellar population (right panel 
in Fig.[8]l. This is due to the effect of the dust attenuation, more 
important for bluer wavelengths where stronger Balmer absorp- 
tion lines are found, especially for young stellar populations. In 
this case, slightly variations in the dust attenuation determination 
can produce the differences in the HafH/3 ratio, more evident for 
small equivalent widths of H/3. We see in Fig.[8]that, in general, 
the relative errors in the determination of the HafH/i ratio are 
similar regardless the dust attenuation fitting, and then, the accu- 
racy in the determination of the gas content is not depending on 
the procedure applied. 

5.2.2. Kinematics 

One of the advantages of using IFU data is the possibility to 
analyse spatially resolved properties of the galaxies. In particu- 
lar, kinematic studies can reveal substructures in galaxies such 
as kinematically decoupled cores in the stellar velocity field and 
twists in the gas velocity field, and also more hidden substruc- 
tures as disks and bars. The instrumental setups used in these 
observations, although limited in spectral resolution, allow some 
studies of this kind. Data obtained with V300 and V600 grism 
allow the measurement of redshifts, while the resolution of the 
V600 grism (FWHM= 5.4 A) will allow measurements of ve- 
locity dispersion above ~ 140 km s"'. In order to asses the 
accuracy in the measurements of redshift and velocity disper- 
sion for these observing setups, we applied the fitting method 
explained above to a set of mock spectra from single stellar pop- 
ulations (an old and a young SSP), and S/N = 50 (per pix) in 
the continuum. In this case, we considered our mock spectra of 
galaxies at different redshifts (from 0.015 to 0.025, in steps of 
0.001) and they were broadened to different velocity dispersions 
(from 150 to 350 km s"', in steps of 25 km s"'). We applied 
the decoupling method in the spectral range AA395Q - 4600 AA, 
which comprises clear absorption features suited for this study. 
In addition, due to the short wavelength considered in this ex- 
periment, a fixed value for the dust attenuation was considered, 
and then, it is not fitted in this case. Finally, for an estimation of 
the random error in these measurements, we simulated a grid of 
spectra of each mock galaxy as follows. An average rms of the 
mock galaxy in areas free of sky lines, and galactic emission, 
or absorption lines was calculated and added by bootstrapping 
method to the SSP-Hgas model spectrum. This was repeated 100 
times, and the fitting was done then on the resulting spectra. The 
results are presented in Fig. |9l where we show the median rel- 
ative errors obtained for the velocity dispersion and the redshift 
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Fig. 9. Estimated relative random errors in the measurements of 
velocity dispersion (left panels) and redshift (right panels) for 
mock galaxies with an old (upper panels) and young (lower pan- 
els) underlying single stellar populations, as explained in the 
text, for observations with the medium resolution grism (V600). 



for each mock galaxy. The error bars represent the standard de- 
viation computed from the simulations. We find that for both un- 
derlying stellar populations, we recover the redshifts with errors 
< 5%. As expected, the relative errors in determination of the 
velocity dispersion are high for values close to the instrumental 
setups. It is possible to measure velocity dispersion with a rela- 
tive error of ~ 30% for cr > 200 km s"' (i.e., intrinsic values of 
the velocity dispersion of galaxies above 140 km s '). 

5.2.3. Ionised Gas 

The Integral Field Spectroscopy allows the study of the resolved 
physical properties of the ionised gas in galaxies, such as (i) 
local nebular reddening estimates based on the Balmer decre- 
ment; (ii) the oxygen abundance distributions based on a suite 
of strong line diagnostics incorporating reddening-corrected line 
ratios (e.g.. Ha, HyS, [On], [Om], [Nii], and [Sii]); and (iii) 
measurements of ionisation structure in H ii regions and diffuse 
ionised gas using the well-known and most updated forbidden- 
line diagnostics in the oxygen and nitrogen lines. For all these 
studies, it is necessary measurements of the fluxes of the emis- 
sion lines in the individual spectra. To test the accuracy in mea- 
sured fluxes for the ionised gas, we followed the procedure ex- 
plained above to create a set of mock galaxies at redshift z - 0.02 
with a velocity dispersion cr = 150 km s"' and S/N = 50 (per 
pix) in the continuum. We measured the main emission lines 
(see § 15.1b in 29 different regions in NGC628, and we created 
29 spectra with gaussians of these values. In this way, the emis- 
sion line ratios of the mock galaxies are realistic and known in 
advance. Then we have applied the decoupling method and mea- 
sured the emission line fluxes. Finally, we have compared these 
results with the input data from NGC628. In Fig. [10] we plot 
the relative random error of the ratio of Hor/HyS as a function of 
the equivalent width for the HjS emission line. We obtain rela- 
tive errors of the ratio of Ha/HyS < 10% when EWh^ > 5 A 
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Fig. 10. Study of the recovering of the properties of the ionised 
gas for the for observations with the medium resolution grism 
(V600). Estimated relative random errors in the measurements 
of the Hcc/HjS ratio as a function of the equivalent width of H/? 
emission lines (EWh/j) on mock galaxies with an old (left pan- 
els) and a young (right panels) underlying stellar population and 
different values for the dust attenuation: magenta for Ay - 0.1, 
red for Al/ = 0.3, green for Ay - 0.6 and blue for Ay = 1.0. Top 
panels are for dust attenuation fixed in the program, and bottom 
panels for dust attenuation fitted in the program. Each block of 
the histogram represents an error of 5%. 



for old stellar populations, while we obtain errors of < 20 % al- 
ways that EWh/3 > 6 a for young stellar populations. There is a 
tendency to over-estimate the emission of H/? for old stellar pop- 
ulations (left panels), while we find the opposite behaviour for 
young stellar populations (right panels). This is expected since 
the Balmer absorption lines are deeper for a young stellar pop- 
ulation and an accurate decoupling is more difficult if the emis- 
sion lines are not very intense. In addition, this effect is clearer 
when the dust attenuation is also fitted in the decoupling method 
(lower panel in Fig.fTOll. 

5.2.4. Stellar Populations 

As explained before, the decoupling method also provides with 
the best combination of stellar populations to reproduce the 
underlying stellar population of galaxies. In this section we 
study the recovering of the luminosity-weighted age and metal- 
licity from the data obtained in this work. With this purpose, 
we applied the decoupling method on spectra of mock galax- 
ies with different underlying stellar populations and gas con- 
tent (see § 15.21 ). at redshift z - 0.02 with a velocity dispersion 
cr = 150 km s-i. We considered different S/N (= 5, 10, 25, 50, 
80, 100 per pixel) in the continuum of the spectra. In addition, 
we created a grid of spectra from each mock galaxies for error 
analysis, following the procedure explained in previous sections. 
Then, we applied the decoupling method in two cases: (i) the 
dust attenuation is a well known and fixed value (results for the 
luminosity-weighted age in Fig. fTTT i. and (ii) the dust attenua- 
tion is a free parameter in the fitting procedure procedure (results 
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Fig. 11. Study of the recovering of the luminosity-weighted age 
as a function of the S/N in the continuum for the V600 grid. In 
this case, the dust attenuation has been considered known and 
fixed when applying the decoupling method. Diff'erent colours 
are used for different values of H/? emission line. 



b) Young SSP 



S/N (pix) 



c) 01% young SSP 



S/N (pix) 



d) 10% young SSP 
I » 



1 



S/N (pix) 



S/N (pix) 



Fig. 12. Same figure than Fig. [TT] when the dust attenuation is 
also computed with the decoupling method. In this case, the es- 
timated random errors (and their standard deviations) are signif- 
icantly higher than in the previous example. 



in Fig. [T2I 1. In both cases, we represent the relative error in the 
luminosity-weighted age for the four cases of underlying stellar 
population. The error bars show the standard deviation for these 
values computed from 100 simulations, as in previous sections. 

As it can be seen, the relative random errors in the 
luminosity-weighted age for mainly old stellar populations (pan- 
els a, c and d in Figs. [TTIandfTSb are ~ 30-55% while the esti- 
mated errors and their standard deviation are higher, in average, 
when considering younger stellar populations (panel b). In ad- 
dition, the relative errors are significantly higher in the case of 
fitting simultaneously the dust attenuation (Fig. \T2\ . It is clear 
from this analysis that the spectral range (similar for data ob- 
tained with V300 and V600 grids) is not wide enough to break 
the age-metallicity-dust degeneracy. Additional information in 
other spectral bands, as ultraviolet and infrared, are desirable in 
order to obtain a reliable value for an independent estimate of the 
dust attenuation. In that case, dust attenuation would be fixed to 
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Fig. 13. Parameters of the most similar SFH spectra (r -left 
panels- and metallicity -right panels) obtained with the decou- 
pling method, as a function of the S/N in the continuum for the 
V600 grid. In this case, the dust attenuation has been considered 
known and fixed in the top panels. Results when the dust attenu- 
ation has been fitted are presented in the bottom panels. Symbols 
and colours are used for different underlying stellar populations 
of the mock galaxies, as it is indicated in the label. 



this value in the fitting process and a more accurate determina- 
tion of the rest of parameters will be obtained. 

In addition, the relative random errors in the luminosity- 
weighted metallicity are extremely high, above the 100% in all 
the cases. This indicates that the grid of templates, based in 
six stellar population mod els, are too simple fo r the estimate 
of metallicity. As noted b y iMacArthur et al.l (l2009h and showed 
bv lSanchez et al.l (l201 Ibh . the analysis of the stellar populations 
is limited by the templates used. Tentative tests with a grid of 
20 and 25 SSPs decrease the error in the luminosity-weighted 
metallicity to ~ 35 %, and therefore a study of the stellar pop- 
ulations is possible with this data by assuming the long compu- 
tational time necessary. It is important to highlight, in any case, 
that the recovering of the kinematics and the properties of the 
ionised gas are well performed with the simple grid considered 
in this paper, and a further analysis of the stellar population pa- 
rameters is feasible with an appropriate grid of templates. 

Due to the uncertainties obtained in the luminosity-weighted 
age and metallicity when a grid of SSPs is considered in the 
decoupling method, we decided to compare our mock galaxies 
with spectra obtained from a given SFH (see § 15. 1.1b . As part of 
the fitting procedure explained in § 15.11 there is an additional op- 
tion to compute the most similar spectrum to the input one from 
a collection of spectra. This option is independent of the global 
fitting to decouple the ionised gas and the stellar population, and 
we can feed the procedure with a different set of template spec- 
tra for this experiment. We used this option to explore which 
SFH spectra (§ 15. 1.1b is more appropriate to describe our mock 
galaxies. In Fig. [T3]we represent the t parameter (left panels) 
and metallicity (right panels) for these mock galaxies, when the 
dust attenuation is fixed (upper panels) of fitted (lower panels). 
The results in both cases are similar, obtaining low values of t 
(4 < T < 6) for predominantly old stellar populations and a high 
value (r ~ 20) for younger stellar populations. In both cases, the 
metallicity is well determined, as expected since we consider 
two metallicity values for this experiment (see § 15. 1.1b . 
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Fig. 14. Comparison of the Ha/H/J, [Niil^eSSS/Ha and 
[Om]/I5007/Hj6 ratios, derived from central spectra for the 
galaxies in this sample and the SDSS spectral, using the multi- 
SSP fitting procedure, plus the comparison between the S/N in 
both spectra. 



5.3. Comparison with the SDSS data 

There are 28 galaxies observed in the sample with available 
SDSS spectroscopy (marked in TablelTJ. In order to compare the 
quality of both data sets, we extracted an integrated spectrum in- 
side a circular aperture of 5" from the datacubes, centred in the 
intensity peak of the galaxy at the V-band. The extracted aper- 
ture is selected to match the aperture of the SDSS fibers, taking 
into account possible errors in the centroid accuracy and see- 
ing effects in the SDSS observations, and the dithering scheme 
adopted in this work. Then we applied the decoupling method 
to the spectra of both samples, deriving their physical proper- 
ties. A full comparison is shown in Fig. [14] Four of the galax- 
ies from SDSS present S/N < 10.0 (per pixel, computed over 
the full spectra), while just one of the galaxies in our sample 
are below this level. As expected, there is a good agreement be- 
tween the values from both samples, especially for the stronger, 
clearly de-blended emission lines, as [Oiii]/15007/HjS (determi- 
nation coefficient for a least squares fits to the data - 0.79). 
For those galaxies with different measurements, it is possible 
that the SDSS fiber would not be well centred in the galaxy, and 
we are just comparing different regions of the galaxy. In addi- 
tion, the S/N of the galaxies observed in this work is higher in 
average than the presented for the SDSS data. 



5.4. Aperture effects in spectroscopy 

One of the advantages of Integral Field Spectroscopy is the pos- 
sibility of integrating different combinat ions of the observed 
spectr a for particular studies. For example. iRosales-Ortega et"an 
(I2OIOI) analysed aperture selected spectra of individual HII 
regio ns withi n the galax ies of the PINGS survey, and 
ISanch ez et all (1201 Ibh (and IViironen et alJl20TTI) explored the 
aperture effects in the galaxies of this survey, using the IFU as a 
large aperture spectrograph. 
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Fig. 16. Same as in Fig. [T4]for the comparison of the values de- 
rived from central and integrated spectra for the galaxies in this 
sample. 



In this section, we analyse two aperture spectra derived for 
each galaxy in the sample: a central spectrum, obtained as ex- 
plained in the previous section, and a 30" diameter aperture 
spectrum, which comprises most of the observed regions of each 
galaxy above the 5cr detection limit. This spectrum is a good rep- 
resentation of the integrated properties of each galaxy. In Fig.fTSl 
we present both aperture spectra for the subset of galaxies in 
Fig- El As it can be seen, the integrated spectra present stronger 
emission lines than the central ones, as expected since they are 
covering also the external parts and different structures of galax- 
ies with usually more amount of ionized gas. There are also dif- 
ferences in the stellar continuum of both apertures for some of 
the galaxies, with redder colours in the central spectra than in 
the integrated one, which point to older stellar populations in 
the central part of the galaxies. For a further comparison, we 
applied the fitting technique presented above to both spectra, de- 
riving the main properties for the galaxies in the sample. We 
represent the values derived from both spectra in Fig.[T6]for the 
main emission line ratios. Again, it can be seen differences be- 
tween both measurements, as expected since they are integrating 
different regions in the galaxies. In particular, the central aper- 
ture is mainly dominated by the bulge of galaxies, where little 
ionised gas is usually found. On the other hand, the integrated 
aperture covers the most of the galaxy size and includes arms, 
bars, and other structures in the galaxies with different gas and 
stellar content. 

As an example of the impact of the differences obtained 
when analysing both spectra, in Fig. [T7] we present the 
classic diagnostic diagram for emission lines galaxies (e.g., 
Veilleux & Osterbrock 1987), using the [Om]/I5007/Hj6 and 
[N ii]/16583/Hq' ratios from central (blu e filled squares) and inte- 
grated (green filled circ les) spectra. T he Kauffmann et al.l (l2003h 
(blue dashed line) and iKewlev et alj (1200 ll) (red solid Une) de- 
marcation curves usually invoked to distinguish between star 
forming and AGN galaxies, respectively, are also plotted. For 
the integrated spectra, 34 of 46 galaxies (i.e., 73.91 % of the 
sample) present accurate emission line ratios (with an error of 
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Fig. 15. Comparison of the central (red) and the integrated (black) spectra for the subset of galaxies presented in Fig.[T] normalized 
to the flux of the central spectrum at 5000 - 6000AA. Differences between both spectra are plotted in green. 



< 30 %). With these data it is found that 30 galaxies are clas- 
sified as pure star forming galaxies (88.23 % of the galaxies 
with accurate emission line ratios), 4 galaxies are in intermedi- 
ate region (11. 76 %) and there are no AGN candidates. 2 galax- 
ies in the sample show little gas (Fua < 0.1 10"'^ cgs Units). 
Based on the properties of the central spectra, just 25 (of 46) 
galaxies present accurate emission lines ratios (i.e., 54.35 % of 
the sample), 14 of them are classified as star forming galax- 
ies (56.00%), 2 galaxies (8.00%) are AGN candidates: IC 2515, 
classifies as Sy2 in the BAT catalogueFI ( Cusumano et al. 201()|), 
and CGCG 428-059, ga laxy already detected in radio and X-rays 
dHardcastle et al.ir2007h . what indicates that probably hosts an 
AGN. There are 9 galaxies in the intermediate region (36.00%). 
In particular, 2MASX J13062093-H5318232, classified as AGN 
in NED database (see Table [1] appears in the intermediate re- 
gion, but the eiTor in the [Nn]/16583/HQ' ratio is > 30 % and it 
will be compatible with its classification as AGN. Finally, in this 
case 18 galaxies show little gas in their central spectra. 

It is clear from this comparison that important differences 
are found when when analysing and comparing spectra extracted 
with different apertures over a galaxy. The two dimensional data 
from IPS surveys can be used to analyse these aperture effects 
(e.g. JViironen et aT1l201 ih . which could be critical for the com- 
parison and interpretation of galaxies at different redshifts. A 
large number of observed objects using this technique will al- 



|http : //www .pa . iasf . cnr . it/cgi-bin/bat/main/fonti . cgi | 



low to obtain empirical calibrations to compare in an appropriate 
way the values obtained in the Local Universe with galaxies at 
higher redshifts. 



5.5. Spatially resolved properties 

The main advantage of using IPS for the study of galaxies 
is the abilit y to derive spatially r esolved spe ctroscopic prop- 
erties (e.g.^ de Zeeuw et al.1 2002; iRo sales-Or tega et alJ I20T0I: 
iBlanc et al.ir2010l: ICastillo-Morales et ahi201Q) . Using IPS, we 
obtain a spectrum for each position of the galaxy, and there- 
fore the decoupling method explained before can be applied for 
each of the individual spectra within the obtained datacubes with 
enough signal-to-noise. In this way, we obtain two-dimensional 
distributions (or maps) of the parameters derived by this analy- 
sis. Although it is out of the scope of this article to analyse in 
detail each of them, we present here a few examples to show the 
quality of the derived maps, focused in this case on the properties 
of the ionised gas. 

Pigure [18] shows the two dimensional distribution of the 
[Nii]/16583/Hq' line ratio derived from the fitting procedures 
described in the text for the sub-set of galaxies in Pig. [T] 
The blue dashed contours shows the continuum intensity at the 
Rc- band. The black solid contours shows the Ha observed 
intensity, starting at 0.3 10"'^ erg s"' cm"^. Notice that the 
data presented in this work enable us to reach sensitivity lim- 
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Fig. 18. Colour map of the [Nii]/16583/Ha line ratio derived from the fitting procedures described in the text, for a sub-set of the 
observed galaxies. The black solid contours shows the Ha observed intensity, starting at 0.3 10"'^ erg s"' cm"^, and with a step 
of 1 10 erg s ' cm"^ for each contour. The blue dashed contours shows the continuum intensity at the Rc- band, starting at 
0.03 10""' erg s"' cm"^, with a step of 0.1 10""' erg s"' cm"^ for each contour. 



its of ~ 10 erg s ' cm ^, fainter than previous imaging sur- 
veys i n the Local Universe (e^g. Perez-Gonzalez et al. 2003j 
Helmboldt et al."2004; Meu rer era02006l:iKennicutt et alj200"8t 
Karachentsev & Kaisin,2010l) . what will allow a more accurate 
determination of the SFR in the Local Universe. In addition, we 
present in Fig. [19] the velocity map of the ionised gas (from the 
Ho- emission), derived with the previously described procedure. 
In most of cases it is seen a clear rotational pattern, mostly due 
to the bias towards late-type galaxies in the presented sub-set. 

Finally, Fig. |20] shows the rotational curves along the 
kinematic major axis, derived from the velocity maps of Ha 
(Fig.[T9]l. The colour code indicates the [Nii]/16583/Ha line ra- 
tio presented in Fig. [18] Values for R25 are taken from the RC3 
catalogue (see Table [TJ- From the sub-sample shown in Fig. [20] 
we can see that it is possible a two dimensional study up to 
~ 1 .0 R25 with the data presented in this work, although in most 
of cases we do not reach so large extension. This easy experi- 
ment shows the importance of the sample selection to accom- 
modate the galaxies in the FOV of the PPAK instrument. Closer 
galaxies with high angular sizes will reach R25 out of the FOV, 
while the S/N of small (and faint) galaxies will be not enough 
for most of two dimensional studies. This work shows, there- 



fore, the importance of the sample selection not only in volume 
and limit magnitude, but also in the apparent sizejjf galaxies as 
previous works showed (e.g.. lBershadv et al.ll2010al) . 

These figures illustrate how suitable is the considered ex- 
perimental setup for the study of these spectroscopic properties 
within the full optical extension of the considered galaxies (char- 
acterised by the continuum emission). A detailed analysis of the 
resolved properties of these galaxies will be presented in forth- 
coming papers. 

6. Summary and Conclusions 

We present in this paper a IFS study of nearby galaxies, exploit- 
ing the capabilities of the wide-field PPAK IFU of the PMAS 
spectrograph mounted on the 3.5 m telescope. In particular, we 
explore the requirements to perform a detailed study of the most 
relevant spatially resolved spectroscopic properties of galax- 
ies in the Local Universe along all their optical extension, fo- 
cused on: (i) The main properties of the ionised gas; (ii) The 
luminosity-weighted age and metallicity of the stellar popula- 
tion; (iii) The main kinematic properties of both the ionised gas 
and stellar population. 



14 



E. Marmol-Queralto et al.: Integral Field Spectroscopy of a sample of nearby galaxies. 






IC2515 Vg^^ 




Ml 








(«> 




. Illy 



A RA (arcsec) 



A RA (arcsec) 



CGCG430-045 V^,,!^ 




— o 





UGC00333 r 




o 


















)^ /A 



































A EA (arcsec) 



A EA (arcsec) 



A RA (arcsec) 









IC3304 T^jj^ 


■ 














o 
o 




















o S 










A DEC 










O 1 


\ 0„ 






















1 . 




1 



A RA (arcsec) 




A EA (arcsec) 



Fig. 19. For the same sub-set of galaxies than in previous figures, colour map of the velocity obtained from the Ha emission, once 
corrected by the corresponding systemic velocity. The black solid contours shows the Ha observed intensity, starting at 0.3 10"'^ 
erg s ' cm"^, and with a step of 1 10 erg s ' cm"^ for each contour, as in Fig. [TS] 



Two sample selections have been explored, based on (1) a 
pure redshift selection and (2) a diameter selected sample, and 
two different instrumental setups has been adopted, covering in 
both cases a similar FOV of ~ 1 arcmin^ and spectral range 
(~ 3700 - 7050 AA), but with different spectral resolutions: 
R ~ 500 and R ~ 900. We show that it is possible to adopt a 
fix instrumental setup and observing strategy to observe a large 
sample of nearby galaxies in a reasonable amount of time (i.e., a 
couple of hours per target, including overheads), covering most 
of the optical extension of these galaxies and wide spectroscopic 
optical range with PPAK. 

We present a detailed description of the data reduction, and 
we demonstrate that the large amount of data produced in this 
kind of observations can be handled if a dedicated automatic re- 
duction pipeline is implemented. In addition, an automatic anal- 
ysis procedure is used to handle the large number of spectra in 
the datacubes derived for each individual galaxy, and to obtain 
the their physical properties. 

We perform a set of simulations to explore the accuracy of 
the physical properties derived from the data. These simulations 
indicate that the adopted observational strategy, and in partic- 
ular, the adopted resolutions and integration times are suitable 
to study both the integrated and spatially resolved spectroscopic 



properties of these galaxies. The simulations show that the spec- 
tral range of the data is not wide enough to determine simultane- 
ously the dust attenuation and the stellar populations, and addi- 
tional information from other spectral bands, as UV and infrared, 
is needed to break the age-metallicity-dust degeneration. The re- 
covering of the properties of the ionised gas and the kinematics 
is well performed with the simple grid of templates used here, 
although the determination of the luminosity-weighted age and, 
especially, metallicity depends on the grid of templates used, as 
expected. 

From the comparison of the galaxy data of this work with 
previously published imaging and spectroscopic data from the 
SDSS survey, we show that it is feasible (1) to obtain spectropho- 
tometric calibrated IFS spectra, and (2) the derived spectroscopic 
properties are recovered equally well than with single fiber spec- 
tra. A showcase analysis of the two dimensional properties for 
a sub-set of galaxies of this sample indicates that it is feasible 
to derive the spectroscopic properties for different apertures, and 
to obtain accurate 2D distributions of these properties within the 
optical extension of these galaxies. 

Finally, we want to highlight that Integral Field Spectroscopy 
constitutes a powerful approach for the understanding of the 
global evolution of nearby galaxies from the analysis of their 
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Fig. 20. Rotational curves derived from the maps shown in the previous figure along the axis of maximum rotation. The colour code 
indicates the [N u]A6583flia ratio, in logarithmic scale. The distances are normalised by the R25 radius. 



resolved physical properties. This study was part of theprepara- 
tory works for a larger IFS survey, the CALIFA surveyQ (~ 600 
galaxies at 0.005 < z < 0.030), which will be an important re- 
source for detailed studies of galaxies in the Local Universe. 
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